Superconducting Gap Anisotropy in Ndi 8 5Ce .i5CuO4: Results from Photoemission 
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We have performed angle resolved photoelectron spectroscopy on the electron doped cuprate su- 
perconductor Ndi.85Ceo.i5Cu04. A comparison of the leading edge midpoints between the super- 
conducting and normal states reveals a small, but finite shift of 1.5-2 meV near the (ir,0) position, 
but no observable shift along the zone diagonal near (tt/2,tt/2). This is interpreted as evidence for 
an anisotropic superconducting gap in the electron doped materials, which is consistent with the 
presence of d-wave superconducting order in this cuprate superconductor. 

PACS numbers: 79.60.Bm, 73.20.Dx, 74.72.-h 



It is now generally agreed that the hole doped (p- 
type) cuprate superconductors have an order parameter 
whose majority component is of d x 2— „2 symmetry. Angle 
resolved photoelectron spectroscopy (ARPES) has pro- 
vided a key piece of early evidence by the demonstration 
of a large momentum space anisotropy of the electronic 
excitation gap pL Quantum interference devices have 
given a definitive attestation to the existence of an or- 
der parameter that changes sign in momentum space Q . 
Corroborating evidence has been provided by a number 
of other techniques Q . 

In contrast, the experimental and theoretical situation 
on the small family of electron-doped (n-type) cuprate 
superconductors Q, such as Ndi.ssCeo.isCuC^, has not 
reached a similar detente. Tunneling microwave pen- 
etration depth H , and Raman j7j] measurements all yield 
evidence for some kind of nearly uniformly gapped Fermi 
surface. In addition, Ndi.gsCeo.isCuC^ is the only mea- 
sured cuprate superconductor that does not exhibit the 
so-called zero-bias Andreev bound state on its (110) sur- 
face that has been interpreted as a consequence of a mo- 
mentum space sign change in the superconducting order 
parameter (|,|§]. 

However, recent work has called this picture of a uni- 
formly gapped Fermi surface into some doubt. There 
is evidence from measurements using scanning SQUID 
microscopy on tri-crystal films, that the electron-doped 
materials have an order parameter with a large d-wave 
component |i"cf| . In addition, there has been specula- 
tion that the ordering at low temperature of the free 4/ 
moments on the rare earth Nd ions affects the usual mi- 
crowave penetration depth measurements by altering the 
low temperature magnetic permeability |@|. A recent re- 
analysis of microwave data that incorporates a correction 
for the Nd spin ordering, as well as experiments on the 
related paramagnetic material Pr2— x 

CezCuC^ have been 
interpreted as being more consistent with a c?-wave sce- 
nario Q. 

On the theoretical side, there have been problems 



reconciling the s-wave experimental results with cur- 
rent models. As in the hole doped case, strong on-site 
Coulomb repulsion in the copper oxygen plane strongly 
favors an anisotropic order parameter that has positive 
and negative lobes. 

We have recently restudied the Ndi.gsCeo.isCuC^ sys- 
tem using high resolution ARPES and found two dis- 
tinct features in the low energy spectra, which contrasts 
with the behavior of the p-types with hole doping, where 
the systematics of the single low energy feature inter- 
polates between the prototypical undoped insulator and 
the optimally doped superconductor JO]. In the case of 
Ndi.85Ceo.i5Cu04, the spectra are dominated by a broad 
highly dispersive feature which resembles that found in 
the undoped insulator Ca2Cu02Cl2 Q. In addition, a 
weaker feature at lower binding energy (< 80meV) is ob- 
served near the (ir, 0) position. These two features should 
not be viewed as deriving from seperate "bands" of inde- 
pendent origin, as they conspire together to give a Fermi 
surface with the expected Luttinger volume. 

In this letter we focus on an energy scale much lower 
than that of our previous work [hjfl (below 5 meV). 
The central intellectual issue is the anisotropy of the su- 
perconducting gap, which should manifest itself in the 
ARPES spectra as a displacement of the leading edge to 
higher binding energy as a gap opens up. 

The single crystals of Ndi.ssCeo.isCuC^ used in this 
study were grown by the traveling solvent floating zone 
method. Details of this growth can be found elsewhere 

s 

ARPES measurements were performed at beamline 5-4 
of the Stanford Synchrotron Radiation Laboratory. This 
system is a Scienta SES 200 electron spectrometer cou- 
pled to a normal incidence monochromator (NIM) that 
is capable of better than 10 meV energy resolution and 
0.2 degree angular resolution. Data was taken in Scienta 
"Angle" mode, in which an entire ±5.5° angular cut can 
be imaged simultaneously. The incident photon beam an- 
gle was at approximately 45 degrees to the sample surface 
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with its polarization in the direction of the angular cuts. 
Samples are positioned by a custom sample manipulator 
that can be cooled below 10 K. The sample tempera- 
ture was measured by a calibrated diode mounted near 
the measurement position. The temperature difference 
between the diode position and the actual measurement 
position was found to be negligible previously from mea- 
surements on another diode that had been mounted to a 
sample holder in a identical fashion to the samples. 

All Ndi.85Ceo.i5Cu04 measurements were taken with 
16.5 eV photons, which has been found empirically to 
give a high photoionization cross-section. Samples were 
aligned by Laue diffraction prior to vacuum chamber in- 
troduction. Shiny flat surfaces resulted from cleaving the 
samples in situ by the top-post method at 10 K. LEED 
was performed after the UPS measurments to check sur- 
face structure and quality. For each measurement set, 
the temperature was cycled a number of times to ensure 
repeatability and an absence of problems from surface 
ageing. No such problems were detected over the typical 
measurement time of 7 hours per k-space cut. To improve 
statistics of these intrinsically low signal experiments the 
same temperature spectra were summed together. 
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FIG. 1. Plot of the LEM around the apparent Fermi surface 
of a few cuprate superconductors of comparable T c 's in units 
of the fractional distance from (tt/2,tv/2) to (tt,0.3it). 

In order to see signatures of superconductivity on a 
system with an expected small gap, it is crucial that 
one compare leading edge midpoints (LEMs) from the 
same sample above and below T c . In the past it has 
been standard practice to simply compare LEMs of the 
spectral weight onset to the Fermi energy of a reference 
gold sample, but this kind of analysis is not reliable in 
the small gap n-type systems. To make this point more 
clear, in Fig. 1 we compare the trend of the LEMs of 
Ndi.85Ceo.i5Cu04 referenced to the Fermi energy to that 
of one plane bismuth based cuprates with similar T c 's as 



our Ndi.85Ceo.i5Cu04 sample. All spectra were taken at 
approximately 11 meV resolution and with the polariza- 
tion parallel to the Cu-0 bond. The bismuth cuprate 
data is adapted from Ref . |lj] . A few distinct trends are 
evident in the hole doped compounds. There is a notable 
decrease in the maximum gap as one proceeds from un- 
derdoping to overdoping. This is consistent with previous 
ARPES results as well as with other techniques. There 
is a definite increase in the LEMs as one goes around 
the Fermi surface in accordance with one's d-wave ex- 
pectation. Similar behavior is seen in YL^CuaOy-z 
and Lai. 85 Sr .i5CuO4. |R||l7j In contrast, this analysis 
which reveals a large anisotropy on the hole doped com- 
pounds, does not reveal any clues as to anisotropy in 
Ndi.85Ceo.i5Cu04 as the LEM referenced to Ep shows 
no systematic trends along its Fermi surface. This is a 
consequence of the fact that the superconducting gap by 
any measure is much smaller than in the hole doped ma- 
terials and hence is obscured due to finite resolution and 
subtle changes in the lincshape around the Fermi surface. 
If the systematics of the lineshape did not change very 
much around , one would still have a good measure of 
relative changes in the gap energy, despite not being able 
to measure its absolute magnitude exactly. This has been 
done, for example, in the imaging of a gap in Bi2212 and 
YBa2Cu307_j; that follows the d-wave functional form, 
despite large measured negative gap energies in the nodal 
region ]l6||lS| ]. In the present case, due to the previously 
mentioned weaker component that appears near (ir, 0) 
the systematics of the lineshape do change and a simple 
comparison with Ep is not valid. One must then compare 
the shift of LEMs from spectra above and below Tc at 
the same momentum space point on the same sample in 
order to accurately measure the gap value. This should 
remove extraneous effects of the final state and otherwise 
eliminate effects (e.g. photoionization matrix elements) 
other than a gap opening up. 

We took coarse energy resolution data at 10K with over 
500 energy distribution curves (EDCs) over a Brillouin 
zone octant and symmetrized across the zone diagonal. 
Fig. 2 shows a 30 meV integration about the Fermi en- 
ergy of these EDCs plotted over a Brillouin zone quad- 
rant. Along the circularly shaped Fermi surface, one sees 
three regions of high intensity. The central one corre- 
sponds to the previously mentioned broad, highly dis- 
persive feature moving into the integration window. The 
other two are the intensity maxima of the smaller low 
binding energy feature. In the inset is a typical low 
energy electron diffraction (LEED) pattern that shows 
bright well formed spots with a symmetry commensu- 
rate with that of the bulk. The extreme surface sensitive 
nature of LEED is an excellent cross-check of the surface 
quality. We took high resolution angle resolved cuts at 
two temperatures through points in two different re- 
gions of momentum space as shown by the arrows in the 
figure ( (0.467T, 0.467r) and (it, 0.277r) exactly, but here- 
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after refered to as near (ir/2,ir/2) and (n, 0.3n) ). In 
both cases the photon energy was 16.5 eV, polarization 
was parallel to the cut direction, and the energy resolu- 
tion was 11 meV. 




30 meV integration about Fermi energy 
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FIG. 2. (color). A 30 meV integration around the Fermi 
energy of each EDC plotted over the Brillouin zone quadrant. 
Two high resolution angle resolved cuts [labeled (a) and (b)] 
were taken as marked by the arrows. These are displayed in 
Fig. 3. In the inset is the LEED pattern of the sample cleaved 
in situ at 10K. 

In Fig. 3a and 3b, we present color scale plots of the 
spectral function in a ±5.5° (±21% n) angular windows 
from the cuts near (ir/2,ir/2) and (ir, 0.3n) respectively. 
In the spectra near (7r/2, 7t/2) we find a dispersion which 
is universal for the cuprates in this region of the Bril- 
louin zone. A large broad feature disperses towards the 
Fermi energy, sharpens, and then disappears as it passes 
above E^. In 3b, which is a cut through kp at approx- 
imately (7r, 0.37r), the spectra are best characterized by 
a large hump feature that disperses only slightly while 
the smaller low energy feature disperses toward E^, loses 
weight, and then disappears. This behavior has been 
characterized in more depth elsewhere (lj]. In Fig. 3c 
and 3d, we compare spectra at 0.5° steps taken at 30K 
(red) and 10K (blue) at a few near kp positions for each 
cut. In 3c from the spectra near (tt/2,tt/2) we see that, 
aside from some small thermal broadening, there is no 
temperature induced change. This is quite different from 
that of spectra from the kp crossing near (7r, 0.37r) where 
there is a systematic displacement by ~1.5-2meV of the 
leading edge to higher binding energy in the supercon- 
ducting state. 

In Fig. 4 there is an expanded image of a spectrum 
with one of the large temperature induced shifts in Fig. 
3d. To more systematically quantify the result, we fit 
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FIG. 3. (color). At the top labeled (a) and (b) are false 
color plots representing the intensity of the spectral function 
in a ±5.5° angular windows (±21% it) near (7r/2,7r/2) and 
(-7T, 0.37r) respectively. In (c) and (d) we compare EDCs taken 
at 30K (red) and 10K (blue) from the fe-E regions represented 
by the horizontal bars in panels (a) and (b) respectively. Ar- 
rows denote the direction that angular cuts are displayed as 
referenced to the arrows in insets of the Brillouin zone. Gap 
values represent the results of the fitting procedures for the 
curves. 

the lineshape with a simple phenomenological model. 
For broad spectral features at the Fermi level, the su- 
perconducting state spectra can be approximated by a 
10K Fermi function multiplied by a linear spectral func- 
tion whose onset edge is displaced from the Fermi en- 
ergy by the superconducting gap energy. Likewise, the 
normal state spectra can be modeled as that of a 30K 
Fermi distribution at finite temperature multiplied by the 
same (but non-displaced) linear spectral function. These 
model spectra are convolved with a Gaussian of FWHM 
of 11 meV to simulate the experimental resolution. The 
agreement between the experimental curves and the fits 
are quite good within this picture with a gap parameter 
of 1.9 meV. 

The above temperature dependent measurements were 
repeated on 7 different samples (4 samples at both kp po- 
sitions, two only at (n, 0.37r) and one just at (7r/2, 7r/2) ). 
In the inset of Fig. 4 we plot the temperature dependent 
shift of the LEM at both kp crossings for all samples. 
Although there is some scatter in the data, one can see 
that, with the exception of one sample, they all show a 
1.5-2 meV shift at the (n, 0.3n) position and a negligible 
one at the (tt/2,w/2) position. This shift anisotropy is 
interpreted as a consequence of the opening of a super- 
conducting gap that is maximal near (n, 0) and minimal 
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FIG. 4. At the bottom of the panel are EDCs from Icf 
near (n, 0.3n) which is curve 2 in from Fig. 3d. They are 
plotted along with a simple fit convolved with the resolution. 
At the top is the same fit unconvolved with the experimental 
resolution. The fitted gap parameter is 1.9 meV. In the inset is 
a plot of the temperature dependence data from all 7 samples 
taken near {-n/2,%/2) and (tt, 0.37r). The data are slightly 
offset from each other in the horizontal direction for display 
purposes. 



termediate temperatures the superconducting gap may 
not have opened fully. In addition, there may be some 
background contribution that partially obscures the gap 
signal. Lastly, it is well documented that photoemission 
analysis based on LEMs consistently underestimates the 
maximum gap value (in some cases by as much as a fac- 
tor of 2) as compared to the intrinsic value defined as the 
quasiparticle energy at the antinodal position. |l§| ]. 

In conclusion, we believe there is evidence for a 
anisotropic superconducting gap in the electron doped 
material that is consistent with a cZ x 2- 2 pairing 
state. Despite the strong differences between the p 
and rt-type compounds in the larger scale electronic 
structure(~leV), their superconductivity appears to 
share the same symmetry and is therefore possibly of 
similar origin. The experimental data was recorded at 
the Stanford Synchrotron Radiation Laboratory which is 
operated by the DOE Office of Basic Energy Science, Di- 
vision of Chemical Sciences and Material Sciences. Addi- 
tional support comes from the Office of Naval Research: 
ONR Grants N00014-95-1-0760/N00014-98-1-0195. The 
crystal growth work was supported in part by Grant-in- 
Aids for Scientific Research from the Ministry of Educa- 
tion, Science, Sports, and Culture, Japan, and the New 
Energy and Industrial Technology Development Organi- 
zation of Japan (NEDO). 



or zero along the zone diagonal. This is consistent with 
the presence of a d m 2- y 2 superconducting order param- 
eter. However, we cannot rule out a very anisotropic s- 
wave order with a small gap along the zone diagonal that 
is below our detection limit (approximately the same as 
the scatter in the data at (ir/2,ir/2), < 0.3?7jeT^). The 
reason for the lack of a temperature dependent shift on 
the one sample is unknown. Its spectral features and 
dispersions were consistent with the others. It is possible 
that it suffered from being in poor thermal contact with 
the cold stage. 

The 2A/ksT c ratio in the narrow doping range of su- 
perconducting n-type materials is typically measured to 
be smaller than those in the hole-doped ones. There is 
no consensus on where they should fit analogously into 
the hole doped phase diagram, but perhaps the T 2 de- 
pendence of their resistivity suggests that they are more 
overdoped-like. We may then expect a low 2A/fcsT c ra- 
tio. The gap value of 1.5-2 meV obtained with our tech- 
nique is consistent with, but slightly smaller than, the 
gap values reported by other techniques There are 

a few possibilities for this. These other measurements 
typically take place at 4.2 K (0.15 T c ), in contrast to 
our measurement at approximately 0.5 T c . At our in- 
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